Condensation Transitions in Nonequilibrium systems 



M. R. Evans 

School of Physics, 
The University of Edinburgh, 
Mayfield Road, 
Edinburgh EH9 3JZ, 
U.K. 

February 2, 2008 

Abstract 

Systems driven out of equilibrium can often exhibit behaviour not seen in systems 
in thermal equilibrium- for example phase transitions in one-dimensional systems. In 
this talk I will review several 'condensation' transitions that occur when a conserved 
quantity is driven through the system. Although the condensation is spatial, i.e. a 
finite fraction of the conserved quantity condenses into a small spatial region, useful 
comparison can be made with usual Bose-Einstein condensation. Amongst some one- 
dimensional examples I will discuss the 'Bus Route Model' where the condensation 
corresponds to the clustering together of buses moving along a bus-route. 

1 Introduction 

Broadly speaking, one can consider two types of nonequilibrium systems: those relaxing 
towards thermal equilibrium and those held far from thermal equilibrium e.g. by the system 
being driven by some external field. In the latter case the steady state of the system will not 
be described by usual Gibbs-Boltzmann statistical weights rather it will be a nonequilibrium 
steady state. A natural way to construct a nonequilibrium steady state is to drive the system 
by forcing a current of some conserved quantity, for example energy or mass, through the 
system. Such systems are known as driven diffusive systems (DDS) j2]. 

In recent years the possibility of phase transitions and phase separation in one- dimensional 
nonequilibrium systems has been explored and some examples are by now well studied. To 
appreciate the significance one should recall the general dictum that in one-dimensional 
equilibrium systems phase ordering and phase transitions do not occur (except in the limit 
of zero-temperature, or with long range interaction) [Hj. 
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Let us briefly review work on one-dimensional phase transitions in driven systems. A 
very simple one-dimensional driven diffusive system is the asymmetric simple exclusion 
process (ASEP). Here particles hop in a preferred direction on a one-dimensional lattice 
with hard-core exclusion (at most one particle can be at any given site). Indicating the 
presence of a particle by a 1 and an empty site (hole) by the dynamics comprises the 
following exchanges at nearest neighbour sites 

10 ^ 01 with rate 1 

01 ^ 10 with rate q (1) 

The open system was studied by Krugjlj and boundary induced phase transitions shown 
to be possible. Specifically one considers a lattice of N sites where at the left boundary site 
(site 1) a particle is introduced with rate a if that site is empty, and at the right boundary 
site (site N) any particle present is removed with rate (3. Thus the dynamical processes at 
the boundaries are 

at site 1 — > 1 with rate a 

at site N 1 -> with rate (3 . (2) 

These boundary conditions force a steady state current of particles J through the system. 
Phase transitions occur when liniAr^oo J exhibits non-analyticities. The steady state of this 
system was solved exactly for the totally asymmetric case E] and more recently for the 
general q case [Tj |S] . When q < 1 the phase diagram comprises three phases: a high-density 
phase where the current is controlled by a low exit rate /3 — one can think of this is queue 
of cars at a traffic light that doesn't let many cars through; a low-density phase where 
the current controlled by a low injection rate a — think of this as a traffic light that does 
not let many cars onto an open road; a maximal- current phase where both a, (3 are high 
(a,/3 > (1 — q)/2) and the current is J = (1 — q)/A. Note that since increasing a and (3 
doesn't increase the current, the current is saturated. In the maximal current phase generic 
long-range correlations exist, an example being the decay of particle density from the left 
boundary to the bulk value 1/2 which is a power law ~ 1/x 1 / 2 where x is distance from the 
left boundary. 

On the line a = (3 < (1 — q)/2 which separates the high and low density phases one 
finds coexistence between a region of low density in the left part of lattice and a region of 
high density on the right separated by a 'shock' where the density changes sharply over a 
microscopic distance. 

Perodic systems (i.e. a ring of sites) can also exhibit phase separation when inhomo- 
geneities or defects are introduced. A very simple example is to introduce into the asym- 
metric exclusion process a 'slow bond' through which particles hop with a reduced rate. 
Then in the steady state one can obtain phase separation between a region of high density 
behind the slow bond and a region of low density in front of the slow bond. Moving defects 
(i.e. particles with dynamics different from that of the others) have also been considered 
and exact solutions obtained fTJElE!- One can think of a slow agricultural vehicle on a 
country road with a large queue of cars behind it and open road in front of it. 
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A further question is whether systems related to the hopping particle models described 
so far, but without inhomogeneities, can exhibit phase ordering. A very simple model was 
introduced in ^3] comprising three species of conserved particles, amongst which all possible 
exchanges are allowed. However a key feature is that the dynamics has a cyclic symmetry 
i.e. A particles move preferentially to the left of B particles which move preferentially to 
the left of C particles which in turn move preferentially to the left of A particles. The 
model exhibits strong phase separation into pure domains of A B C. Similar strong phase 
separation occurs in other related models j!5j . 

A final class of transitions in one-dimensional hopping particle models is that involving 
spatial condensation, whereby a finite fraction of the particles condenses onto the same 
site. Examples include the appearance of a large aggregate in models of aggregation and 
fragmentation ^H] and the emergence of a single flock in dynamical models of flocking |17j . 
We will analyse a simple example of a condensation transition which occurs in the zero-range 
process which we now define. 

2 The zero-range process 

The zero-range process was introduced by Spitzer [0] and recent applications and develop- 
ments have been reviewed in [10]. We consider a one-dimensional lattice of M sites with 
sites labelled /i = 1 . . . M and periodic boundary conditions (more generally one can con- 
sider the zero-range process on a lattice of arbitrary dimension). Each site can hold an 
integer number of indistinguishable particles. The configuration of the system is specified 
by the occupation numbers of each site fi. The total number of particles is denoted by L 
and is conserved under the dynamics. The dynamics of the system is given by the rates at 
which a particle leaves a site \x (one can think of it as the topmost particle — see Figure la) 
and moves to the left nearest neighbour site \x— 1. The hopping rates u(n) are a function of 
n the number of particles at the site of departure. Some particular cases are: if u(n) = n 
then the dynamics of each particle is independent of the others; if u(n) = const for n > 
then the rate at which a particle leaves a site is unaffected by the number of particles at 
the site (as long as it is greater than zero). 

The important attribute of the zero-range process is that it has a steady state described 
by a product measure. By this it is meant that the steady state probability Pdn^}) of 
finding the system in configuration {rii, n 2 . . . um} is given by a product of factors /(n M ) 



Here the normalisation Z(M, L) is introduced so that the sum of the probabilities for all 
configurations, with the correct number of particles L, is one. 




M 



(3) 
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In the basic model described above, f(n) is given by 

n 1 

f(n) — 1 I — r fc> r n > 1 
y ' 11 u(m) 

m=l v ' 

= 1 for n = (4) 



To prove ()3I4|) one simply considers the stationarity condition on the probability of a 
configuration (probability current out of the configuration due to hops is equal to probability 
current into the configuration due to hops): 

^ 9{n^)u{n li )P(ni . . .n^, . . .n L ) =) j 0(n^)u{n^ + i + l)P(ni . . . n M -l, n M+i +l . . . n L ) ■ 

A* A 1 

(5) 

The Heaviside function highlights that it is the sites with n > 1 that allow exit from 

the configuration (lhs of but also allow entry to the configuration (rhs of ©)• Equating 
the terms // on both sides of and cancelling common factors assuming (jHJ), results in 

«W/( n M-i)/( n A*) = + !)/( n M ~ ^/('V+i + !) ( 6 ) 

This equality can be recast as 

u ( n v) jf~^ L T\ = u ( n v+± + 1)^-77-^—^— = constant (7) 
Setting the constant equal to unity implies 

and iterating (jHJ) leads to (jlj). 

We can easily generalise to consider an inhomogeneous system by which we mean the 
hopping rates are site dependent: the hopping rate out of site \i when it contains n M particles 
is u^in^). It is easy to check that the steady state is simply modified to 

1 L 

p (K}) = ^X)II^W ( 9 ) 



where /» are given by 



U(n) = I f — p—r for 71 > 1 

= 1 for n = (10) 



The proof is identical to that for the homogeneous case, with the replacement of ^(n^) by 
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Figure 1: Equivalence of zero range process and asymmetric exclusion process. 



There exists an exact mapping from a zero-range process to an asymmetric exclusion 
process. This is illustrated in Figure 1. The idea is to consider the particles of the zero-range 
process as the holes (empty sites) of the exclusion process. Then the sites of the zero-range 
process become the moving particles of the exclusion process. Note that in the exclusion 
process we have M particles hopping on a lattice of M + L sites. A hopping rate in the zero 
range process u(m) which is dependent on m corresponds to a hopping rate in the exclusion 
process which depends on the gap to the particle in front. So the particles can feel each 
other's presence and one can have a long-range interaction. 



3 Condensation Transitions 

We now proceed to analyse the steady states of form (jHJ) and the condensation transition 
that may occur. The important quantity to consider is the normalisation Z(M, L) as it 
plays the role of the partition sum. The normalisation is defined through the condition 

M 

z(m,l)= (11) 

n\,ri2...riM M 

where the 8 function enforces the constraint of L particles. The normalisation may be 
considered as the analogue of a canonical partition function of a thermodynamic system. 
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We define the 'speed' v as the average hopping rate out of a site 



1 n 

V = Z(M L) ^ 5 (}2 n ^- L ) u ( n i)\\f^ n ^) 

Z(M,L-1) 

~ Z(M,L) [ ' 

where we have used fl9H0|) . Note that (fT^j) tells us that the speed is independent of site 
and thus may be considered a conserved quantity in the steady state of the system. In 
the totally asymmetric system considered in Section |2] the speed is equal to the current of 
particles flowing between neighbouring sites. The speed is a ratio of partition functions of 
different system sizes (II 2|) and corresponds to a fugacity. 

We now use the integral representation of the delta function to write the partition func- 
tion as 

Z(M,L)=f£:z-W IIW. (13) 

where 

oo 

F,{z) = * m Ufa) . (14) 

m=0 

For large M, L ffTTTf) is dominated by the saddle point of the integral and the value of z at 
the saddle point is the fugacity. The equation for the saddle point reduces to 

T M A 

M M ^ dz M ; y J 

which, defining <fi = L/M, can be written as 

* = 5^W' (16) 

In the thermodynamic limit, 

M -> oo with L = (pM , (17) 

where the density is held fixed, the question is whether a valid saddle point value of z 
can be found from (|Tfi]). We expect that for low (j) the saddle point is valid but, as we shall 
discuss, there exists a maximum value of z and if at this maximum value the rhs of (|T6|) 
is finite, then for large <ft (fTK|) cannot be satisfied. We now consider how condensation may 
occur in the inhomogeneous and the homogeneous case. 
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3.1 Inhomogeneous case 



To give an idea of how a condensation transition may occur we consider the case u M (m) = 
for m > i.e. the hopping rate does not depend on the number of particles at a site. / M is 
given by 

' 1 



and the probability of occupancies {rii, n 2 , . . . ,Um} is 



1 M ( 1 

p «"^ "" |)= «ny • a*) 

The mapping to an ideal Bose gas is evident: the L particles of the zero-range process are 
viewed as Bosons which may reside in M states with energies determined by the site 
hopping rates: exp(— fiE^) = Thus the ground state corresponds to the site with the 

lowest hopping rate. The normalisation Z(M, L) is equivalent to the canonical partition 
function of the Bose gas. We can sum the geometric series ()14j) to obtain F^ and F^ then 
taking the large M limit allows the sum over /i to be written as an integral 

duV(u) (20) 



where V(u) is the probability distribution of site hopping rates with w min the lowest possible 
site hopping rate. Interpreting V{u) as a density of states, equation (|2Uj) corresponds to the 
condition that in the grand canonical ensemble of an ideal Bose gas the number of Bosons 
per state is <fi. The theory of Bose condensation tells us that when certain conditions on the 
density of low energy states pertain we can have a condensation transition. Then (jlfij) can 
no longer be satisfied and we have a condensation of particles into the ground state, which 
is here the site with the slowest hopping rate. 

A very simple example is to have just one 'slow site' i.e. u\ = p while the other M — 1 
sites have hopping rates u^—l when fi > 1. Using the mapping to an exclusion process, 
this corresponds one slow particle i.e. agricultural vehicle example described earlier. One 
can show JT] that for a high density of particles in the zero range process (low density of 
particles in the corresponding asymmetric exclusion process) we have a condensate since 
site 1 contains a finite fraction of the particles. In the low density phase the particles are 
evenly spread between all sites. 



3.2 Homogeneous case 

We now consider the homogeneous zero-range process where the hopping rates u(n) are site 
independent. Then (fTljl is independent of \i and reads 

oo n 

*w = e n 

n=0 m=l 



mm) 



(21) 
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The fugacity z must be chosen so that F converges or else we could not have performed 
(114)1 . Therefore z is restricted to z < (3 where we define (3 to be the radius of convergence 
of F(z). From (|2*T|) we see that (3 is the limiting value of the u(m) i.e. the limiting value of 
the hopping rate out of a site for a large number of particles at a site. We interpret (fTH|) as 
giving a relation between the density of holes (number of holes per site) and the fugacity z. 
The saddle point condition (JT5)l becomes 

, zF'(z) , , 

*=TW < 22) 

Given that the rhs of ()2*2*)) is a monotonically increasing function of z we deduce that density 
of particle increases with fugacity. However if at z = (3, the maximum allowed value of z, 
the rhs of (|22|) is still finite then one can no longer solve for the density and one must have 
a condensation transition. Physically, the condensation would correspond to a spontaneous 
symmetry breaking where one of the sites is spontaneously selected to hold a finite fraction 
of the particles. 



Thus, for condensation to occur (i.e. when </> is large enough for (j22j) not to have a 
solution for the allowed values of z) we require 

lim ^-fl < oo . (23) 
F(z) v ; 

We now assume that u(n) decreases uniformly to (3 in the large n limit as 

u(n) = (3(l + ((n)) (24) 
where ((n) is a monotonically decreasing function. Analysis of the series 

OO ( 71 

F((3) = J>xp j - ^ In [1 + C(m)] 

n=0 L m=l 

co ( n 

F'{(3) = 5>expi-^ln[l + C(m)] \ (25) 

n=0 I m=l J 

reveals that the condition for condensation is simply that F'(f3) is finite and this occurs if 
u{n) decays to f3 more slowly than /3(l + 2/n). (This is easiest to see by expanding In [1 + (] 
and approximating the sum over m by an integral in ()23|) .) 

It is interesting to translate this result into the language of the exclusion process. In this 
context we can have condensation if the hop rate of a particle into a gap of size n decays as 
(3(1 + 2/n) therefore there is an effective long range interaction. 



3.3 Bus route model 

As an example of this let us consider the 'bus route model' [IB]. The model is defined on 
a Id lattice. Each site (bus-stop) is either empty, contains a bus (a conserved particle) or 
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contains a passenger (non-conserved quantity). The dynamical processes are that passengers 
arrives at an empty site with rate A; a bus moves forward to the next stop with rate 1 if 
that stop is empty; if the next stop contains passengers the bus moves forward with rate (3 
and removes the passengers. 

The model thus defined has not been solved but simulations reveal two regimes. At high 
bus density the gaps between buses are evenly distributed. However at low bus density 
there is a condensed regime where the lead bus has a large gap to the next bus in front of 
it with bus-stops full of passengers in between. The other buses have small gaps between 
them. Thus the buses form a jam of buses and after a long delay all arrive at a bus-stop at 
once. 

The bus route model can be related to the zero-range process by a mean-field approx- 
imation in which we integrate out the non-conserved quantity (passengers). The idea is 
that a bus-stop, next to bus 1 say, will last have been visited by a bus (bus 2) a mean 
time ago of n/v where n is the distance between bus 2 and bus 1 and v is the steady state 
speed. Therefore the mean-field probability that the site next to bus 1 is not occupied by 
a passenger is exp(— Xn/v). From this probability an effective hopping rate for a bus into a 
gap of size n is obtained by averaging the two possible hop rates 1, (3: 



We can now see that this mean-field approximation to the bus-route model is equivalent to 
a homogeneous zero-range process discussed earlier. 

Since u(n) decays exponentially the condition for a strict phase transition in the thermo- 
dynamic limit is not met. However on any finite system for A sufficiently small, an apparent 
condensation will be seen. In the bus route problem this corresponds to the universally 
irritating situation of all the buses on the route arriving at once. 

4 Conclusion 

We have shown how the zero range process exhibits two kinds of condensation transition. 
One is due to having an inhomogeneous system i.e. we get condensation of particles onto 
the site with the slowest hopping rate. Although the condensation is spatial the mechanism 
is equivalent to Bose condensation in an Ideal Bose Gas. 

The other type of condensation occurs on a homogeneous systems and involves the 
spontaneous selection of a site onto which a finite fraction of the particles condense. Recently 
this condensation mechanism has been used to understand the existence or non-existence of 
phase separation in a general class of one dimensional driven systems |19| I20j. 




(26) 
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